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Abstract. In this paper we analyze urban spatial segregation phenomenon in terms of the 
income distribution over a population, and inflationary parameter weighting the evolution of 
housing prices. For this, we develop a discrete, spatially extended model based in a multi— agent 
approach. In our model, the mobility of socioeconomic agents is driven only by the housing 
prices. Agents exchange location in order to fit their status to the cost of their housing. On the 
other hand, the price of a particular house changes depends on the status of its tenant, and on 
the neighborhood mean lodging cost, weighted by a control parameter. The agent's dynamics 
converges to a spatially organized configuration, whose regularity we measured by using an 
entropy-like indicator. With this simple model we found a nontrivial dependence of segregation 
on both, the initial inequality of the socioeconomic agents and the inflationary parameter. In 
this way we supply an explanatory model for the segregation-inequality thesis putted forward 
by Douglas Massey. 

Mathematics Subject Classification: 91D10, 91B72 



1. Introduction 

The spatial structure of a city is the result of a wide and complex set of factors. The way in 
which different economic activities and social groups spread over the urban space is the matter of 
different and complementary theories in Sociology, Geography, Politics and Economy |5ll9"l ll5IIH4| . 
Housing patterns can be understood as the result of the complex interrelationship between individ- 
uals' actions constrained by social, political and economical rules [HlESj- Residential segregation 
is the degree to which two or more groups live separately from one to another in different parts 
of the urban space p. 282]. This phenomenon is concurrent to several social problems as the 
concentration of low opportunities to get a well earned job, low scholar development in children 
of segregated areas, premature parenthood between young people, and the emergence of criminal- 
ity [HIES- 

Residential segregation has been the subject of extensive research in social sciences for many years. 
It is a multifactorial phenomenon mainly determined by socioeconomic factors like race and income 
distribution, as well as factors associated to the structure of the urban space 1131 1171 1211 12"2"] . 

There are two main quantitative approaches to segregation, the phenomenological one, which 
relies on segregation measures and indexes ^01 D2I E3 EH|, an d the theoretical one, based on 
computational or mathematical models [31 [HI ESI CZI CHI EI3 EU E21 EDI EH E21 EH] 
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In this paper we propose a spatially extended model based in ideas of Portugali 28 and Schelling 30 
1311 13*2"] , to study the relationship between income inequality and residential segregation. The basic 
thesis proposed by Massey and colleagues [2J is that the degree of spatial segregation experi- 
enced by a society increases with its level of inequality. This relationship has been reformulate by 
Morrison |24j as the segregation-inequality curve. 

The organization of the article is as follows. In the next section we briefly expose the main 
theoretical studies concerning the residential segregation, emphasizing only the ideas and concepts 
that are relevant to our research. In section 3, we present the mathematical model and the tools 
needed to analyze our numerical simulation, which we do in section 5. Previous to this, in section 
4 we study of asymptotic behavior of the model, and derive some theoretical estimates which we 
consider for the numerical study. Finally, we conclude with a discussion about the potential of our 
model as an explanation tool. 

2. The segregated city 

Residential segregation is a complex phenomenon with several dimensions of analysis, whose gov- 
erning mechanisms are hard to identify. In a first approximation we can however assume that the 
phenomenon is governed by a set of structural and behavioral rules which determine the possibility 
of one individual to get a particular kind of house in a specific location of the city. Since those 
rules are no evident, simplifying hypotheses are required. 

One point of view, based in human ecology, postulates that residential segregation occurs because 
individuals in a city are in mutual competition for the space and its resources. According to 
this approach, competition is the main force driving the residential segregation p. 86]. The 
outcome of this competition is determined by the ability of individuals to struggle for advantageous 
locations in the urban space, i.e. their dominant capacity, which is constrained by sociocultural 
and socioeconomic rules pp. 85-88]. 

There are three main hypothesis about the sociocultural rules governing the residential segregation. 
The first one concerns with the class-selective emigration from poor regions. In a region where 
coexist both poor and less poor people, the latter tend to emigrate to a more wealthy region. This 
mechanism tends to isolate and concentrate poor people, increasing in this way the poverty rate 
of the region. The second hypothesis establishes that neighborhood concentration of poor people 
reflects the general poverty of the urban area. When the average shows a downward trend, neigh- 
borhood poverty rates increase. Finally, the third is related the the racial segregation experienced 
by poor people. Racial bias causes racial segmentation of the urban housing markets, which con- 
curs with high rates of poverty in specific ethnic groups to concentrate poverty geographically |221 
pp. 426-428, and references therein]. These hypotheses are complementary, and were developed to 
explain segregation in north-american cities, where they have been tested. Perhaps in the Latin- 
american case racial and sociocultural factors have a less relevant role, making possible to build 
an explicative model over socioeconomic considerations only. 

Taking into account that "markets are not mere meetings between producers and consumers, whose 
relations are ordered by the interpersonal laws of supply and demand" p. 1], we can formulate 
socioeconomic rules as market mechanisms. The housing market is formed by two kinds of agents: 
residents which are interested in the social and individual value or use of the land commodity, 
and the entrepreneurs which are interested in the exchange value of the land. There is a natural 
conflict between these two of valuations of land. 

There is a set of structural factors that are relevant to housing market dynamics: a) the house- 
building industry; b) the government's housing policy; c) the structure of the property of land; 
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Figure 1 . Segregation-inequality curve according to Morrison et al. HU 



d) the actual spatial structure of city, i.e. location of labor area, residential areas, and trade- 
commerce areas; and e) the income structure of the society. The last one has been considered as 
the most significant for the residential segregation phenomenon. Indeed, Urban Economic theory 
explains the formation of segregated cities through two main arguments. The first one establishes 
that a population of households with heterogeneous income competing for the occupation of urban 
land traditionally results in an income-based stratification of the urban space according to the 
distance to the city center pp . The second one links the concentration of low incomes households 
in some areas, to the existence of local externalities like ethnicity. As a consequence of this, there 
is a households' preference to live in relative homogeneous neighborhoods with respect to either 
income or ethnic similarity [51121132]. 

Two levels of analysis may be considered in the study of residential segregation. At the macro 
level, several structural transformations in the society (changes in income level, tendency to racial 
exclusion, levels of social integration, etc.) are assumed to determine the spatial concentration of 
poverty. This is the level of analysis in [U EH E2 . At the 

micro level, specific discriminatory 
individual behaviors related individual characteristics (sex, age, religion, ethnic group, nationality, 
etc.) influence the choice of a place to live. This is the point of view in EEl EH E01 03 E21 ) and 
it is also the one we adopt here. 

Our model was developed with the purpose of studying Massey's thesis [21 P- 400], which re- 
lates the degree of spatial segregation experienced by a society, to its inequality degree (income 
disparity). This thesis was reformulated as the segregation-inequality curve (see Figure^) by Mor- 
rison 24 j. More precisely, we intend to determine the relationship between these two quantifiable 
phenomena, inequality and segregation, in a situation where the whole dynamics is governed basic 
rules of socioeconomic nature. 
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3. Model Structure 

Our model is inspired in the works of Schelling and Portugali et. al [2S| . Like Portugali's 

models, the physical infrastructure of our "simplified city" is modeled by a two-dimensional lattice 
of finite size A m n = {1, 2, . . . , m} x {1, 2, . . . , n}. A two-dimensional integer vector x = (xi, X2) £ 
h~m,n represent spatial coordinates. At each time step, the price of the house at location x is a 
positive real number. Each house is occupied by a householders or agent, who can be distinguished 
only by his/her socioeconomic status. We quantify this status with a real number taking one of 
three possible values p < m < r (which stands for poor, middle class, and rich respectively). At 
time t, each agent occupies one specific house of the city. Houses are identical in their characteristics 
but differentiable by their prices. 

There are two main mechanism setting up the house's prices dynamics: the neighborhood influence, 
and the householder's economic status. On the other hand, the agents change position subject to 
availability, under the pressure of their housing situation. Agents can move inside the city to 
achieve an optimal match between their status and the price of the house they inhabit. Agents try 
to live in houses with prices, but not too much, that their economic status. 

At time t, the prices of houses are encoded inamxn matrix V 1 := (V*(x)) xe A m n , while the spatial 
distribution of agents is stored mamxn matrix A 1 := (A*(x)) xg A m „ with values in {p, c, r}. The 
price of that house at time t + 1 is given by: 

n ) v t+1 M = A t+1 M + a £|x(i)- y (i)i+|x(2)- y (2)i<2 v \y) 

(> (> 1 ' #{yeA m ,„: |x(l)-y(l)| + |x(2)-y(2)|<2} 

The parameter A weights the influence of the mean price on the neighborhood over the house price, 
and can be thought as an inflationary parameter: the larger A the higher the asymptotic mean 
value of the houses in the city. 

The distribution of the agents at time t + 1 differs from its distribution at time t only by a single 
place exchange between two agents of different status, i. e., at time t we randomly choose two sites 
x, y € A mt n and define 

(2) 5(x,y) := (A*(x) - ^(x)) 2 - (A*(y) - ^(x)) 2 + (A*(y) - ^(y)) 2 - (A'(x) - ^(y)) 2 . 

This quantity should be interpreted as the economic improvement due to the house exchange 
between agents a locations x and y. At each unite time, only two agents can participate in such a 
house exchange. We have A t+1 (z) = A t (z) for z ^ {x, y} and 



(3) 



A t+1 (x) = I 



A*{y) if<5(x,y)>0, 
A*(x) otherwise, 



and similarly for y. The economic improvement due to a house exchange leads to the reduction 
on the economic tension 



(4) T(A,V) = \A-V\:= J £ (A(x) - ^(x)) 2 

x£A„ 



Agents try to minimize economic tension generated by difference between housing price and status. 
In our model this difference plays the same role as the dissatisfaction or the unhappiness in the 
Schelling's model H3IU]. 
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3.1. Analytical Tools. The segregation-inequality curve gives the relation between two charac- 
teristics of the system: inequality and segregation. The income distribution can be interpreted 
as a probability vector. The proportions of the population in each income group would be the 
probability for an individual to belong to that income group. With this idea, income inequality 
can be measure by using Theil's inequality index |351 pp. 91-96]. Adapting this index to our 
situation, we define inequality ofamxn agents' distribution A as follows: 



where for i — p,m,r, q% = #{x £ A m) „ : A(x) — i}/(m x n) is the proportion of agents in 
each income group. The inequality so defined is an entropy-like indicator taking values in the 
interval [0, log(3)]. The minimum I corresponds to the case where the total population is equally 
distributed among all the income groups, and maximum to the limit case of a single income group 
concentrating the whole population. This last limit would be obtained from distributions where a 
given income group includes most of the population. 

The other characteristic we need to determine the segregation curve is the spatial segregation itself. 
Though measures of segregation have been the subject of several works in sociology ^0 p. 283 
and references therein] , it is more suitable to our approach to quantify this characteristic by using 
the degree of order of a given spatial distribution. The idea is to associate the maximum degree 
of order to a spatial distribution which can be easily described, like a single cluster or a periodic 
distribution. On the other hand, a random distribution would have a low degree of order. We 
use an image segmentation technique, the bi-orthogonal decomposition, to associate a degree of 
disorder (the entropy of the bi-orthogonal decomposition), to a given spatial distribution which we 
treat as an image joj p. 131]. Thus, to quantify the segregation (ordering) of the city, we compare 
the entropy of the bi-orthogonal decomposition of the ordered distribution with that corresponding 
to a random distribution. 

The bi-orthogonal decomposition is the bi-dimensional generalization of the Karhunen-Loeve 
transform, but with the advantage that it is sensible to changes in the spatial structure of one 
image. To a bi-orthogonal decomposition it is associate an entropy, which measures the amount 
of information of the image [SJ p. 133]. The interpretation is that in an orderer image, i. e. an in- 
homogeneous distribution of pixels showing a spatial pattern, has low entropy, while to a random 
distribution of pixels corresponds highest entropy 

Consider a m x n positive matrix U, which is supposed to codify an image or a two-dimensional 
distribution. The associate covariance matrix Q := Wll is symmetric, and hence its eigenvalues 
{Xi : i = 1, . . . , m} are real, and the corresponding eigenvectors {fa : i = 1, . . . , m} define an 
orthonorma! basis. The bi-orthogonal decomposition allows us to rewrite the matrix U as, 



where ipi = Ufa for i = 1, 2, . . . , m. The contribution of the submatrix ipifal to the sum is of the 
order of the corresponding eigenvalue |A$|. The information of the matrix U is concentrate in the 
submatrices associated to eigenvalues with the highest absolute value. For positive U, we neglect 
the largest eigenvalue Ai, which can be associated to an spatially homogeneous mode. Using the 
eigenvalue structure of the bi-orthogonal decomposition, we define the information contents of U 



(5) 



1(A) = log(3) + 




in 



(6) 




by 



(7) 
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where 

(8) Pi = J^U- 

,Lfe=2 \ 

For an agent distribution A, the segregation index is 

(9) Sbo(A) = E(H BO ) - H BO {A), 

where E(iJeo) is the expected value of H B o with respect to a random distribution of agents in 
the m x n-dimensional lattice. We have numerically found that E(i?Bo) ~ l°g(3/5 x min(m, n)) 
for min(m,n) < 1000. 

4. Asymptotic behavior of the model 

In order to understand the asymptotic behavior of the model, let us rewrite Q and© in matrix 
form as follows: 

(10) V t+1 = A' + XDV* 

(11) A t+1 = P*A* 

Where both V* and A* have to be considered as m x n-dimensional vectors, D is the averaging 
matrix whose action is defined by , and P is a permutation matrix which permutes at most two 
coordinates. If for the chosen coordinates x,y G A m „ we have 5(x, y) > 0, then P* is the matrix 
permuting those coordinates, otherwise P is the identity matrix. 

After a sufficiently large number of iterations, name T, the agents achieve a spatial distribution 
A* which cannot be improved. From that point on, the distribution of housing prices follows the 
affine evolution V T+t = A* + XDV T+t ~ l , so that 

(12) V T+t = {XD) l V T + (^2(^D) s ^j A* = (XDYV T + (id - (AD)*) (Id - AD)" 1 A*. 

The long-term distribution of housing prices is therefore V* := (Id — XD) 1 A*, so that the eco- 
nomic tension associate to the asymptotic distribution of agents is 

(13) T*(A*) := T(A*,V*) = \A* - (Id - XDy 1 A* \ = A x | (Id - XD)' 1 (DA*)\. 

Because of the non-deterministic nature of the evolution of our system, the asymptotic distribution 
A* is not uniquely determined by initial conditions. Nevertheless, it has to satisfy the following 
"variational principle" 

(14) T*(A*) = nun (p - (Id - AD) -1 ) A* 

where the minimum is taken over the set of all two-sites permutations. This is equivalent to say 
that asymptotically no location exchange can diminish the economic tension. 

4.1. Critical A. For A small, a spatially disordered initial distribution of agents ^4° remains un- 
changed, and the system evolves following an affine law, converging to A* = A°,V* = (Id — 
AD) -1 )^4°. For each spatially disordered initial distribution A , there exists a critical value A c for 
which A a evolves towards an spatially organized state A* . This distribution is composed by rela- 
tively small number of clusters, each one of them consisting of a nucleus of rich agents surrounded 
by middle class ones, while the lower class agents occupy the space left by the clusters. In our 
numerical experiments, which we describe below, we have found that A c essentially depends on the 
proportion of r, m, and p in A . 
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A two sites permutation A* i— ► PA* produces a change in the economic tension 

T*(A*) >-> |A (Id - XD)~ 1 (DA*) + (A* - PA*)\. 

According to Equation (|14() . this change does not make the economic tension to decrease. In order 
for this to be so, it is necessary that 

(15) \PA* - A*\ 2 > 2A((Id - XD)- 1 (DA),PA* - A*). 

for each two sites permutation A* i— ► PA* . We may decompose DA* = A* 1 + /, as the sum of a 
constant vector and a fluctuating one. Since 

(M-XD)-'1 = T ^ J 

and (1,PA* — A*) = 0, then Equation i|T5|) can be written as 

\PA* - A*\ 2 > 2A((Id - XD)- 1 /, PA* - A*). 
Taking this into account, we may define A c as, 

l\ 1 A ((Id - XD)~ 1 f, PA* - A* 

(16) A c = min < A > : - < max ■ u 



v 2 p \PA* - A*\ 2 

where the maximum is taken over the two-sites permutations. For P interchange coordinates x, y, 
we have, 

A ((Id - AP)- 1 /, PA* - A*) _ ((Id - AD)- 1 /)x - ((W - XD)- 1 f) y 



(17) 



\PA* - A*\ 2 2(A*-A*) 



4.2. An a priori estimate for A c . A reasonably good estimate for A c can be obtained as fol- 
lows. Considering A* as a random field, the Central Limit Theorem ensures that, with very high 
probability ZAA* £ [A* - 2a, A* + 2a] , where 



Here N is the number of sites in the computation of the local mean DA*., which in our case is 25. 
Hence, with very high probability, 



max|/ x -/ y | = 4a = - J Pr (r - A*) 2 + Pm (m - A*) 2 + Pp {p - A*) 2 , 

where p r , p m , p P are the proportions of r, m and p in A* respectively. By using the upper bound 

max Xjy |/ x - / y | 



max|((Id-A J D)- i /)x-((Id-AP)- i /) y | < .. N 



we obtain 



A ((Id - XD)- 1 !, PA* — A*) X 2 J Pr (r - A*) 2 + p m (m - A*) 2 + Pp {p - A* f 
max — : / — ^— < x 1 



p \PA* - A*\ 2 ~1-A 5(m-p) 

According to Eq. I|lt)|l . for A > A c we have 

1 A 2 y/p r (r- ~A*) 2 + Pm (m- A*~) 2 + p p (p - A*~) 2 



< z r x 



2 1 - A 5 (m - p) 

Taking this into account, we propose 

(18) A*- 5{m - p) 



4 ^p r (r - A*) 2 + p m (m - A*) 2 + p p (p - A*) 2 + 5 (m - p) 



cS 
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as an estimate for A c . 

4.3. An priori estimate for Sbo(A*) at A c . Suppose that A* is composed by N c clusters of 
comparable size, and suppose also that all of them are symmetric with respect to the coordinate 
axes. In this case A* has a cluster decomposition 

(19) A* :=gQ4, 

i=l 

where QiP} corresponds to the i-th cluster. The vectors Pi and Qi have a belled form with 
maximal value at coordinates where the cluster is located. The vectors {Pi}^ 1 span a vector 
space of dimension n c < N c , for which {pk}kLi ^ s an orthonormal bases obtained from {Pj}^ by 
the Gram-Schmidt process. Using this orthonormal base, we can rewrite A* := Y^kLi 1kp\. , where 
the vectors qk are obtained form Qi after the change of basis. With respect to this new basis, A* 
can be considered a random n c dimensional matrix, therefore Hbo(A*) ~ log(2 * n c /3), and hence 

(20) S B o(A*) » log(min(m, n)/n c ). 

Let ^V m i n be the cardinality of the less numerous class of agents, i. e., 

N mhl = min #{x <E A TOi „ : A*(x) = i}. 

i—p.ra^r 

For small values of the inequality index, i. e., for large N m { ni we have numerically found that 
n c ~ y/N m i n . This is consistent with an agents' distribution A* formed by N c oc -/V m j n clusters, each 
cluster containing nearly the same number of agents of the less numerous class. The corresponding 
cluster decomposition would be obtained from a collection {Qi}™^ of n\ oc %/AWn bell shaped 
vectors, and another collection {Pj}"=i of ri2 ~ V N m - ln nearly orthonormal bell shaped vectors, as 

ni ri2 

In this case we have N c = n\ x ri2 cx N m i n , and n c = n-i ~ v^min- Hence, for low values of the 
inequality index /, we may expect 

(21) Sbo04*) « S* := log(min(m, n)/ y/N~). 

Note that iV m i n does not determine the value of /. If this value is large enough, the variability of 
iV m j n inside the collection of agents' distributions with the same / value produces a large dispersion 
in S* . For this reason it is not possible to define S* as a function of /, therefore there is not a 
unique inequality-segregation curve. 

5. Numerical Results 

We performed a set of numerical experiments in the n X n lattice, for n — 64 and 128. Each 
lattice node represents a house location in our virtual city, where agents and values of houses are 
distributed. At time t, theses distributions are codified bynxn real valued matrices. The agent's 
matrix A* takes only three values, r = 1, m — 1/2 and p = 1/10, representing the income of 
rich, middle class and poor agents respectively. The distribution of house prices V* is a positive 
matrix, which at time t = takes values in the interval [0, 1]. Both spatial distributions, V* and 
A*, evolve interrelatedly according to Eqs. and 0. The unevenness of an agent's distribution 
A is measured by using Theil's index 1(A) defined in Eq. Since this indicator depends only 
on the proportions of rich, middle class and poor agents, then I(A t ) — I (A ) for all t > 0. 
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Modeling the demographic composition of a city in a developing country, we have considered 
"demographic scenarios", i. e., a given number N p of poor, N m of middle class, and N r of rich 
agents N r , such that N r < N m < N p . In order to obtain an inequality-segregation curve, we have 
chosen a one parameter family of demographic scenarios where N p increases from half to the total 
population, while the ratio a = N r /N rn < 1 is kept constant. In this way we obtain the family of 
demographic scenarios 

T a :={(N p ,N m ,N r ) ~ (7/n 2 ,(l-a)(l-ry)n 2 ,a(l-r/)n 2 ) : 1/2 < 1}, 

which we call "regular" . In our numerical simulation we found no qualitative difference among 
different regular families. Instead, we have contrasted the behavior shown by !F a , with that of 
the family of demographic scenarios obtained as follows: for each value of / we choose, among the 
demographic scenarios with this inequality index, the most probable one subject to the condition 
N p < N m < N r . For this we first consider the uniform distribution in {(N p , N m , N r ) : 1 < N < 
n 3 , i — p,m,r}, then under conditions J2i= P m r — n 3 , N p < N m < N r , and I(A*) — I, we 
distinguish the most probable one with respect to the conditional probability distribution. In this 
way we obtain a family f mp of demographic scenarios which we call "most probable" . 

For the regular family and for each lattice size n — 64 and 128, we have considered 10 different 
demographic scenarios. The demographic scenario thus determines the value of the inequality 
index I. These demographic scenarios have been chosen in order to have 10 nearly equally spaced 
values for I in the interval of possible values [I a , log(3)] . Note how the minimum possible value 
I a := log(3) + l/21og(l/2) + a/21og(a/2) + (1 - a)/21og((l - a)/2), depends on the parameter 
a of the regular family. Now, for each demographic scenario we have chosen 5 different values for 
the parameter A around our estimate A*. Once fixed the demographic scenario and the value of 
A, we performed 20 experiments in both, the 64 x 64 and the 128 x 128-lattices. The purpose of 
these experiments was: 1) to determine A c and compared it to our estimate; and 2) to compute 
the value of the segregation index at A = A c . By doing so, we were able to determine a functional 
relation between the inequality index and the segregation index, which allows us to draw the 
segregation-inequality curve. The same experimental protocol was implemented for the family of 
most probable scenarios. In this case we produce 9 demographic scenarios, corresponding to 9 
equally spaced values of I in the interval [0, log(3)]. 

Each experiment started with spatial distributions V° and A randomly generated. The entries of 
V° were always taken independent and uniformly distributed in the interval [0, 1], while for A the 
agents determining a given demographic scenario were randomly distributed in the lattice. The 
experiment consisted in the iteration of the evolution rule, Eqs. Q and (J3J), until a stationary 
distribution A* was reached. In Figure we show the asymptotic distributions of agents A*, ob- 
tained from the same initial condition A by using 4 different values of A. Note how the asymptotic 
distribution acquires a more regular structure as we increase A. 

In order to determine A c , for a given initial agents' distribution A , we compute the evolution of the 
bi-orthogonal decomposition entropy H^o{A t ) 1 considering increasing values of A. If A is small, 
H-Qo(A t ) remains practically constant along the evolution, while for sufficiently large values of A, 
this entropy undergoes a monotonous decreasing until a definite time that we call segregation time, 
at which it attains its limiting value. We illustrate this in Figure where we show H-qo^A*) for 
A < A c and A > A c . Hence, the segregation time may be considered infinite for small values of A, 
and taking finite values for A > A c . We determine A c corresponding to a initial agents' distribution 
^4°, by computing £/bo(^4') for A = A* — 0.05, A*, . . . , A* + 0.15, and taking the smallest of these 
values for which segregation time is smaller than the empirically determined convergence time. 
The time T an initial configuration ^4° needs to attain the asymptotic distribution A*, increases 
with both I(A°) and A. Nevertheless, for the 64 x 64-lattice and A < 1, this time never exceeded 



10 



ANTONIO AGUILERA AND EDGARDO UGALDE 




-1 = 0.8 A = 0.9 



Figure 2. Asymptotic agent's distribution for different values of A in the 64 x 64- 
lattice. Here I (A*) » 0.25 and A c « 0.6 



1200 iteration in our simulations. For the 128 x 128-lattice, this convergence time at A c was always 
smaller than 2000, and never exceeded 3200 for the other values of A. The convergence time T 
appears to increase in proportion to nlog(n), where n is the lattice size. 

In Figure^ we show the behavior of A c as the inequality index changes, for the regular family in 
the 64 x 64-lattice, and we compare it to the behavior of our a priori estimate A*. In Figure |5] we 
display the same comparison for to the 128 x 128-lattice. Since A* depends only on the proportions 
of rich, middle class and poor agents, it is lattice size independent. According to our numerical 
results, our a priori estimate is a reasonably tight lower bound for A c . Let us remark that the 
numerical value of A c slightly depends on the initial condition of the experiment, hence we plot 
the mean value of A c and the corresponding error bars. 

For the family f mp , we show in figures El and \7\ the behavior of A c as a function of the inequality 
index, in the 64 x 64 and 128 x 128-lattices respectively. We compare this to the behavior of 
our a priori estimate A*. Since A c slightly depends on the initial condition, we show A c with the 
corresponding error bars. Once again, A* is a reasonably tight lower bound for the actual value of 
Ac 

In figures 0D and we plot the segregation index Sbo as function of the inequality index, for 
the regular family in the 64 x 64 and 128 x 128-lattices respectively. These are the inequality- 
segregation curves our model produces. We compare those curves to our a priori upper bound 
estimate S* :— log(n/y/Nr) with n = 64 and 128, which we derived at the end of Section 4. In 
figures HHIandlrTlwe plot the same data, corresponding to the family of most probable demographic 
scenarios. Our numerical results show that, in the case of the regular family, our prediction holds 
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Figure 3. Evolution of the segregation index 7Jbo(^')- We show two exper- 
iments in the 64 x 64-lattice, for an agents' distribution with inequality index 
I = 0.5 and A* « 0.68. The green horizontal line corresponds to A = 0.65, the 
blue line was computed with A = 0.75, and the red was obtained using A = 0.8. 




i 



Figure 4. In red A c with its error bars, also as function of Theil's inequality 
index /. In blue A* also as a function of /. Both curves correspond to the regular 
family, with a = 0.4, in the 64 x 64-lattice. 



for inequality indices in the interval < I < 0.5. For the family of the most probable demographic 
scenarios, the a priori upper bound holds up to I = 0.7. The exact value of 5bo(^*) depends on 
the initial condition A , therefore we show its mean value with the corresponding error bars. 
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Figure 5. In red A c , with its error bars, as function of Theil's index I. In blue 
A* also as a function of /. Both curves correspond to the the regular family, with 
a = 0.4, in the 128 x 128-lattice. 



0.95- 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

I 



Figure 6. In red A c with its error bars, as a function of /. In blue A* also as a 
function of Theil's inequality index i\ Both curves correspond to the the family 
f m „ in the 64 x 64-lattice. 



6. Discussion and Conclusions 

The model proposed here is intended as a tool towards the explanation of the well-known and 
accepted thesis of Massey about the relationship between income inequality and spatial distribution 
of individuals in a city. Build upon basic rules governing the dynamics of house prices and house 
holders, this model is able to produce spatial patterns whose degree of order grows with the 
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Figure 7. In red A c with its error bars, as a function of Theil's inequality index 
/. In blue A* also as a function of /. Both curves correspond to the family J- mp 
in the 128 x 128-lattice. 

3r 
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Figure 8. In red the segregation index Sbo with its error bars, as function of /. 
In blue S* = log(64/-\//Vr) also as a function of /. These inequality-segregation 
curves correspond the T a family, with a = 0.40, in the 64 x 64-lattice. 

income inequality of the virtual city. Thus we have a mathematical model of segregation, showing 
a behavior in accordance to the Massey's thesis. This model introduces a control parameter which 
we relate to the adjustments of the house prices during the evolution, in a way that segregation 
occurs only for large enough values of this parameter. 

The model does not use a sophisticated housing prices theory, but simple interaction rules similar to 
those proposed by Portugali and Benenson. Regardless to its validity inside the economic theory, 
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Figure 9. In red the inequality-segregation curve, Sbo as function of/. Here the 
segregation index Sbo is plotted with its error bars. In blue S* — log(128/-\//Vr) 
also as a function of I. Both curves correspond to the T a family, with a = 0.40, 
in the 128 x 128-lattice. 
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Figure 10. In red the segregation index Sbo, with the respective error bars, as 
function of /. In blue S* — log(64/-\//Vr) as a function of Thcil's inequality index 
/. Both curves correspond the J- mp family in the64 x 64-lattice. 



we believe this is a good first approximation for the housing prices dynamics. The mechanism 
responsible for the segregation is a location exchange process similar the one proposed by Schelling. 
This process is driven by the economic tension resulting from the difference between the "economic 
power" of the house holders and prices of the houses they occupy. 
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Figure 11. In red the segregation index Sbo, with the respective error bars, 
as function of Theil's inequality index /. In blue S* = log(128/-\//Vr) also as a 
function /. Both curves correspond to the jF mp family in the 128 x 128-lattice. 

Though our model retains many of the Schelling's ideas, the main difference is nature of the location 
exchange decision. In Schelling's original model, this decision is taken based on a satisfaction 
function that takes in account the number of neighbors of the same kind around a given agent. In 
our case, an agent decides to exchange is location according to the difference between its income 
and the price of its house. We think that this is a more realistic situation in the sense that that 
economic tension can be measure in real life. This mechanism is also better suited for a model of 
a society where ethnic differences are less determinant than differences imposed by the income. 

From the theoretical point of view, the model is interesting by its own. It exhibits an order-disorder 
phase transition via clustering produced from a very simple exchange mechanism. It also appears 
to present some finite size scaling behavior which would be pertinent to further explore. 

Concerning the indicators we use, Theil's is a widely accepted inequality index, with the advantage 
that it is easy to compute and has a direct interpretation. Our segregation index, on the other hand, 
have never been used in this context. Besides this indicator, we previously tried other entropy-like 
measures of the degree of order in a spatial distribution, as well as some direct clustering measures. 
We chose the entropy of the bi-orthogonal decomposition because it is almost as easy to compute 
as the Theil's index, and it has a natural interpretation as the information contents of a picture. 
We are convinced that other segregation measures commonly used in the sociological literature are 
no suited for our purposes, mainly because they are built from a previous organization of the data, 
e. g., the census tracks. In our case, the structuration of the urban space is a result only of the 
location exchange dynamics. 

Summarizing, our model is built from a very simple location exchange rule, based on interaction 
between agents through an economic tension produced by the difference between the agents' eco- 
nomic capacity and the prices of the houses they occupy. This simple mechanism is sufficient to 
produce spatial segregation in accordance to Massey's thesis: segregation is an increasing function 
of the inequality. With this we were able to furnish key components of an explanatory model of 
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spatial segregation in a situation where the economic factors are more significant than either ethnic 
or cultural ones. 
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